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This article gives a formula that allows accurate values of the cumulative 
noncentral chi-square distribution to be computed. Although this distribution 
has been recognized jor a long time, none of the standard references give 
formulae that are suitable for computing accurate values over an extensive 
range of the parameters] approximations in terms of the chi-square dis- 
tribution are usually recommended. A program written by the author, based 
on the formula given here, has been successful for computations involving 
more than 10,000 degrees of freedom. Since many steps are required when the 
degrees of freedom are as large as this, the program is not "fast" but it is 
believed to be accurate. 

I. INTRODUCTION 

The Non-Central Chi-Square Distribution is encountered in many 
statistical problems, one of the most important in communications 
studies being the detection of signals in noise using a square-law 
detector. 1 Marcum discussed this application but concluded that a 
satisfactory algorithm for computing system performance could not 
be based on the formula he used. 2 This article shows that a satisfac- 
tory algorithm can be based on the formula that Marcum derived if 
the expression is expanded in a power series and the terms are properly 
grouped before being evaluated. 

More recently Urkowitz 3 discussed detection system performance 
in which the above distribution arose and recommended that approxi- 
mations in terms of the chi-square distribution, given by Patnaik, 1 
be used for computation. While these approximations are adequate 
for some purposes, it is desirable to have a reliable and accurate 
method of computing values, if only to check the approximations. 



* This work was supported by the U. S. Navy under contract N00039-68-C-3584. 
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H. INTEGRATION OF NONCENTRAL CHI-SQUARE DISTRIBUTION 

If the signal-to-noise power ratio is x for the sum of -q independent 
samples of the output of a square-law detector, the following integral* 
gives the probability that the sum will be y or more. The variables 
are normalized to the variance of the individual noise samples, so 
the average signal-to-noise power ratio for one sample is xf-q and the 
average output per sample is y/tj. Considering one sample of noise 
to be the sum of the squares of two independent gaussian variables 
of unit variance, the integral is related to the noncentral chi-square 
distribution by the conversions given in equation (8). 

Q = | [-) exp (-z - aO/^pCa)*] dz. (1) 

From (Ref. 4, Section 8.445) J_i [2 (zx)*} is the modified Bessel 
function 



I f t )=f flflr" . (2 ) 



Thus 
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(3) 



Notice that 



/: 



exp (-2) 2 "- 1 <fe = r(77, y) (4) 

the incomplete gamma function (Ref. 5, Section 6.5.3) . Since 

f f exp (-0 dt = -f exp (-0 + p f° f" 1 exp (-*) d« (5) 
equation (3) can be written 



+ ^{r(^y) + (^)y'- I ex P (-y) 



♦Sometimes called the generalized Marcum Q-function. See Ref. 2. 
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Summing the terms by columns gives 

_ r(n, y) ir 1 exp (-y) \y , ) f* £ 

Q ~ r« + r<*) L e p ( } & H 

^f + if + ^ ^Sn 
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+ • • ■ and so on (7) 

A satisfactory computing algorithm can be based on equation (7) 
where we notice that Q can be expressed as the sum of two parts, Qi 
= r(7/, y)/T{ v ) which is independent of x, and another part which 
we call Q 2 . 

III. DISCUSSION 

The noncentral chi-square cumulative distribution can be written 
Q(x 2 \v, A) (see Ref. 5, Section 26.4.25), where the distribution is in- 
tegrated from x' 2 to infinity, the number of degrees of freedom is v, 
and the noncentral parameter is A. This integral is the same as that 
given in equation (7) if we put 

v = 2t\ 

X' 2 = 2y (8) 

A = 2x 
so that 

Q(2y | 277, 2x) = Q 1 + Q 2 (9) 

= Q(2y | 2,) + Q 2 
where Q(2y\2 v ) = Q(x 2 |"), the cumulative chi-square distribution 
(see Ref. 4, Section 26.4.2) . 
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If M independent samples of the output of a square-law detector 
are averaged, when the input is narrowband gaussian noise plus a CW 
signal at the center of the band, Q can be used to find the probability 
that a threshold value will be exceeded. Expressing all parameters in 
units of the narrowband noise power, the desired threshold is y/rj, xf-q 
is the signal-to-noise power ratio, and M = 77. 

It is interesting that the Rayleigh distribution, and the Rice dis- 
tribution, are equivalent to the chi-square and non-central chi-square 
distributions respectively, when the latter are expressed in terms of a 
parameter x equal to the square root of x z , and 9 =s 1. 

Marcum 2 gave an expression of the form shown in equation (3) 
for the output of a square-law detector. He stated that it could only 
be used satisfactorily for values of 77 up to about 10. More recently 
Urkowitz 3 has discussed the integration of a square-law detector out- 
put and recommends that the noncentral chi-square distribution be 
computed using an approximation given by Patnaik 1 in terms of the 
chi-square distribution. Patnaik compares with exact values some 
results computed using the approximation and finds errors of the 
order of 1% around Q = 0.5. The accuracy is much less for values 
around unity and for values less than 0.01. 

Brennan and Reed have shown that, when the order of the Bessel 
function in equation (1) is zero, corresponding to one sample, a 
straightforward recursive method applied to the resulting equation 
(6) can be used to compute the integral. They suggested that a 
similar procedure could be used even on the form of equation (1) 
given here. However, as pointed out by Marcum, such a technique 
rapidly becomes useless as 77 increases above about 10. 

A program written by the author, based on equation (7), has been 
used satisfactorily for 77 as large as 8192, and simultaneously for 
values of z/77 up to 0.1. The exact values given by Patnaik were 
checked. Further checks were made possible by the development of a 
uniform asymptotic expansion by S. O. Rice, with which it is possible 
to get results outside the useful range of the algorithm given here. 7 
Table 1 compares values obtained with the author's program 
(CHISQ) and corresponding values supplied by S. O. Rice using his 
uniform asymptotic expansion (UAE), with results obtained using 
the Patnaik 1 and Gauss approximations (Ref. 5, Section 26.4.29). 

The accuracy of the algorithm given in equation (7) decreases as 
.T./77 increases in the table, and the value in the last entry depended 
quite sensitively on the last digit of a 18 digit double precision con- 
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stant used in the program. Notice that even for the last entry, the 
value actually computed, (CHISQ) , appears to be correct to 14 places 
after the decimal point. 

IV. EXTENSION TO A MORE GENERAL INTEGRAL 

A more general integral is obtained by writing, for example, the 
/?th moment of the partial noncentral chi-square distribution, 

Q t = /"**(fj " exp (-z - x)^-i[2(zx) J ] dz. (11) 

The corresponding form of equation (4) is 

I" z exp (-zy- 1 dz = T(S,y) (12) 

where 

i = v + ft (13) 

and the corresponding form of equation (7) becomes 

Q = ' £ ^exp(-x) 1 F 1 (ii; v ;x) 

y*- 1 exp (-y) Vy f* x_ (g + l) r -, 

VKV + 1) r= 2 ri (y + 2) r _, 



+ £ exp(-x) ± *- itt Sd 

I" , ( , + 1)(l? + 2 ) ex P * ^ r 4j r! (, + 3) r _, 



■]■ 



+ • • • and so on ■ (14) 

The confluent hypergeometric function iFi(a; b; x) (Ref. 5, Section 
13.1.10) is denned by 

j? (n . h . „\ _ i i *EL _i_ «(a + l)s a , a(a + !)(« + 2)s 3 
,i ^ ' ' ; " + 61! + 6(6 + 1)2! + 6(6 + 1)(6 + 2)3! + " ' 

Equation (15) conveniently gives an example of Pochhammer's sym- 
bol (a) r (Ref. 5, Section 6.1.22), also used in equation (14). 



CHI-SQUARE DISTRIBUTION 207 

The structure of equation (14) is closely related to that of equa- 
tion (7), so it can form the basis for a useful algorithm to compute 
the integral given in equation (11) . 
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